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ABSTRACT 

We develop a pscudo power spectrum technique for measuring the lensing power spec- 
trum from weak lensing surveys in both the full sky and flat sky limits. The power spec- 
trum approaches have a number of advantages over the traditional correlation function 
approach. We test the pseudo spectrum method by using numerical simulations with 
square-shape boundary that include masked regions with complex configuration due 
to bright stars and saturated spikes. Even when 25% of total area of the survey is 
masked, the method recovers the E'-mode power spectrum at a sub-percent precision 
over a wide range of multipoles 100 < ^ < 10''. The systematic error is smaller than 
the statistical errors expected for a 2000 square degree survey. The residual i?-mode 
spectrum is well suppressed in the amplitudes at less than a percent level relative to 
the E'-mode. We also find that the correlated errors of binned power spectra caused 
by the survey geometry effects are not significant. Our method is applicable to the 
current and upcoming wide-field lensing surveys. 
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1 INTRODUCTION 

Large scale structure deflects light rays as they prop- 
agate from distant galaxies to us, thus distorting the 
shap e s of these galaxies ( e.g., iBartelmann fc Schneiden 
I2OOII : iHoekstra fc JainI l2008l . for thorough reviews). This 
weak lensing or cosmic shear signals measures a combi- 
nation of the total matter distribution projected along 
the line-of-sight and the angular diameter distance. 
Since the first measurements of weak lensing only a 



decade ago (IVan Waerbeke et al.l 


2OOOI: Wittman et al.ll2000l: 


Bacon. Refregier fc EUid I2OO0I: 


Kaiser, Wilson fc LuDDind 


2000f). there have been maior improvements as surveys 



continue to grow in size and depth (e.g. iFu et al. 





2008 




201c 



llchiki. Takada fc Takahashil |2009| : ISchrabback et al 
for the latest results). Gravitational lensing is one of the 
most promising methods of constr aining cosmology includ - 
ing the nature of dark energy (e.g. iTakada fc Bridl3l2007h . 
There are various on-going and planned surveys aimed at 
studying dark energy through the high-precision weak lens- 
ing measurements: the CFHT Legac y Surve^l^, the Hype r 
Suprime-Cam Weak Lensing Survey ( Mivazaki et al.|[2006l ). 



the Dark Energy Survey (DES a^nd ultimately Lar ge Syn- 
optic Survey Telescope (LSSTfl Euclid ( Rcfrcgi cr et al.l 
l2010t ). and Joint Dark Energy Mission (JDEM). 

How should we analyze these new lensing data set? 
Most researchers use the two-point correlation function to 
characterize the cosmic shear signals. The correlation func- 
tion method can be easily applied to complex survey ge- 
ometries involving partial sky coverage and masked regions. 
However, the errors in the measur ement are highly cor- 
related between different bins (s ee ISchneider et al.l |2002| : 
Ijoachimi. Schneider fc Efleill200a . for the detailed studies). 
Even if the shear field follows Gaussian statistics, which is 
a good approximation in the linear regime, there are large 
correlations between different angular-scale bins. These cor- 
relations are even larger, and more model-dependent on the 
smaller scales that contain most of the current observational 
information. Since the accurate estimate of the covariances 
is essential for robust cosmological constrai nts, a large num- 
ber o f numerical simulations are necessary (|Semboloni et al] 
l200'if ). 

The power spectrum, the Fourier- or Harmonic- 
transformed counterpart of the two-point correlation func- 



^ http:/ /www. cfht.hawaii.edu/Science/CFHLS/ 



^ http:/ /www. darkenergysurvey.org 
^ http://www.lsst.org 
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tion, is an alternative means of measuring the cosmic shear 
correlations. Whilst the correlation function and power spec- 
trum are mathematically equivalent, the power spectrum 
measurement of cosmi c shear has been used less (see the 
COMBO 17-survey bv lBrown 6^11120031 '). The power spec- 
trum approach has a number of advantages: its theoretical 
interpretation is simpler and there are weaker correlations 
between band powers at different multipoles. For example, 
the different bins are independent for the Gaussian field or 
on large angular scales. Even for small angular scales af- 
fected by nonlinear structure formation, the power spectrum 
covariances are relatively well understood through both ana- 
lytic al models and simu l ations of nonlinear structure forma- 
tion (IHu fc Whitell200ll:ICoorav fc Hull200ll : iTakada fc JainI 
I2OO9I : ISato et al.ll20Q9l : iPielorz et al.ll2009l ). The disadvan- 
tage is the presence of finite sky coverage and masked re- 
gions, which breaks the orthogonality of Fourier/Harmonic 
components. One needs to properly deal with the survey 
geometry effect to estimate unbiased power spectrum. 

The purpose of this paper is to eliminate this dis- 
advantage. We employ the pseudo power spectrum tech- 
nique, which is well developed in the CMB studies (e.g. 
IWandeh. Hivon fc Gorskill200ll ) Fl. For the first time, we ap- 
ply the method to recover the lensing power spectrum from 
the shear field taking into account incomplete survey geom- 
etry. To assess the performance of this method, we make 
simulated shear maps including a realistic configuration of 
masked regions due to bright stars and saturated spikes. Fur- 
thermore, we develop the method for both the full-sky and 
flat-sky approaches. The full-sky approach is adequate for 
reconstructing large angular-scale modes that are relevant 
for the curvature of the sky. On the other hand, the flat-sky 
approach should serve as a practically useful approximation 
of sub-degree scale modes, which carry most of useful cos- 
mological information in the shear power spectrum. We find 
that the pseudo power spectrum method allows for an un- 
biased estimate of the underlying i5-mode power spectrum 
over a range of angular scales we study. We also show that 
the residual B-mode power spectrum, which is leaked from 
i?-mode power due to an imperfect reconstruction, can be 
well suppressed. Our method can be applied to the existing 
data and forthcoming weak lensing surveys. 

The paper is organized as follows: Section [2] describes 
the pseudo spectrum method to deconvolve shear power 
spectra with inhomogeneous survey mask. Section [3] de- 
scribes the simulation maps we use to test the deconvolution 
method. We employ two different simulation maps: one is 
Gaussian shear maps and the other is the ray-tracing simu- 
lations of shear maps including the non-Gaussian effects due 
to nonlinear structure formation. SectionUshows the results 
of both the full-sky and flat-sky approaches. Section [S] is de- 
voted to the summary and conclusions. 



2 METHODOLOGY: RECONSTRUCTION OF 
SHEAR POWER SPECTRUM 

In this section we briefly review a method for reconstructing 
shear power spectra from the pseudo-spectrum estimators. 



We take into account an imperfect survey geometry due to 
survey boundary and masking effect. The method is analo- 
gous to the one used in e stimatin g CMB pol arization power 
spectra (|Kogut et al. I lioba : .Brown. Castro fc Tavlor.,2005i V 



2.1 Full Sky Formalism 

Since the shear field is a spin-2 field, the E- and B-mode 
harmonic coefHcients of the shear fields 7i(i = 1,2) can be 
expressed in the spherical harmonic expansion as 



Elm ± iBi„ 



dfifi [7i(n) ± ij'zin)] ±2y4(n) 



and the inverse relation is 

71 (n) ± 172(6) = ^ [Ei,n, ± iBim] ±2yim(n) 



(1) 



(2) 



where ±2Yim. is the spin-2 spherical harmonics, and n de- 
notes the unit vector specifying the angular direction on the 
sky. The integration range is over the full sky. 

In the linear regime, the statistical information in the 
map is fully encoded in the power spectra 



ct 
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cf 
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(3) 
(4) 
(5) 



In the single lens limit, the shear field arising from a scalar 
gravitational field should be a gradient or curl-free field 
{Bim = 0). The multiple lensing effect generates B-mode 
power spectrum, but their power is ~ 10* times smaller than 
the E-mode power. Thus the i?-mode power spectrum effec- 
tively contains all information on the cosmic shear, i.e. E- 
mode power spectrum is equivalent to that of the projected 
mass density along the line-of-sight between source galaxies 
and an observer. Hence the B-mode can be used as a moni- 
tor of residual systematic effects. The standard methods to 
separate E/B mode correlation functions involve integrals of 
the measured correlation functions down to arbitrary small 
scale or up to very large scale. As the scale range accessi- 
ble from a fini te sky data is limited, residual uncertainties 
are generated (Schneider ct al." 199^: 'Critte nden et al.ll2002l . 
also sec Schneider, Eiflcr & Krausc 2010 for a new method 
using the limited-range integration of correlation function 
to separate the i5-mode). 

Observational effects, such as a finite sky coverage and 
bright star masks, limit the survey area to a region K{n). 
The observed shear field is modified as 



71 (n) ± 172 (n) = _K'(n)(7i(n) ± 272(6)). 



(6) 



Without weighting, K{h) = if the position vector fi lies 
in masked regions or regions outside the survey, otherwise 
K{n) = 1 within the survey. This finite sky coverage couples 
modes and generate artificial B modes. We can describe the 
observed shear fields in terms of "pseudo E and B modes", 
denoted as Eim and Bim 



* See ISeliakI l ll998l') and iHu fc Whitl hOOTh for the maximum 
likelihood method of shear power spectrum estimation. 



Elm ± iBu 



df^fi [7^(n)(7i(n) ± ^72(6))] ±2YiUn). (7) 
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These pseudo E and B modes are related to the true E and 
B modes as 



II' mm' ■) 



(8) 



through a convolution kernel, 



±2 =f2 



47r 

/ I' I 



II 



m m m 



,(9) 



where ( ] are known as the Wiener 3 7 sym- 
\^ mi 7712 ms J 2, 

bols (see the references in lDahlen fc Tromdll99^ ) and Kim 

is the harmonic transform of the mask function K{h): 



Kim = i dni,K{h)Yi*m{n)- 



(10) 



The pseudo power spectra Cf^, C^^ and Cf'^ are de- 
fined similarly to the equations ©-((S]) using pseudo E and B 
modes (eq. [8]). After straightforward algebraic calculation, 
one can find that the pseudo and true spectra are related to 
each other via 



Ci=J2^ii'Fi'Ci'+'Ni, 



(11) 



where we have introduced the vector notations C; = 
{Cf^jCP^jCf^) and so on for notational simplicity, and 
Fi is the pixel window function. In the above equation we 
include the shot noise contribution arising from the intrinsic 
ellipticities of source galaxies. The intrinsic noise is simply 
modeled as the convolved noise power spectrum N; with 
the pixel window and the mask. Non-zero components of 
the mode-mode coupling matrix M;;' are given as 



21' + 1 

Stt 



^(2r-hl)K;r'[l + (-l)'+''^ 



I I' I" 
2-2 



(12) 



ji.fEE,BB _ j,tBB,EE 



IV 

21' + 1 
Stt 



^(2r 4- l);Cr'[l - (-1)'+''+'"] 



I I' I" 

2-2 



(13) 



M ,- 



3,EB _ 21' + 1 



Y^{21" + l)K,v 



I I' I" 
2-2 



with ICi being defined as 



21 



l-^'^KimKl 



(14) 



(15) 



Note that M^^, 



= M,7, Equation dTT]) tells 



equivalently a leakage of i?-mode power into the B-mode 
even if Cf ^ = 0. 

The underlying power spectra can be reconstructed by 
solving the equation (|lip inversely. The resolution in multi- 
pole space is lim ited by the survey area, i.e. the finer binning 
less than l{ = y^Tv/faky does not improve the statistical sig- 
nificance of the power spectrum reconstruction. Since the 
lensing power spectrum does not have fine scale structures 
in multipole space, such a coarse binning is sufficient to cap- 
ture the shape of shear power spectrum. Also the coarse bin- 
ning significantly reduces the computational cost of measur- 
ing the power spectrum in a wide range of multipoles. We 
therefore measure the binned power spectra defined as 



Cb = y ^ PbiCi, 



(16) 



where the index "fe" denotes the &-th multipole bin, and 
'Y^t'' represents the summation over I between l'"^^^ and 
^miV — 1. Here we use a binned operator Pbi so that the 
binned power becomes the average of dimensionless power 



over I between I'^^l^ to l''^^'^ — 1 



,(6) 

+ 1 



'min Snin 

The deconvolved binned spectrum is obtained as 

leb' 

Cb = {m-^)bb' Pb'iiCi - (n,)mc). 

The mode mixing matrix for the binned spectra is 

teb I'eb' 

Mbb' ^Y^^Yl ^u'F^Qi'b'^ 
I I' 

where Qw is the reciprocal of Pm 



(17) 



(18) 



(19) 



(20) 



i{i + iy 

The equation (jlSp is the key equation for reconstructing the 
shear power spectra from masked shear maps in the full-sky 
approach in Section ID 

2.2 Flat-Sky Approximation 

In this subsection, we present a formalism of a pseudo- 
spectrum method in the fiat-sky approximation. The fiat-sky 
approximation is enough applicable to current lensing sur- 
vey, such as the CFHT survey covering a sky of about 2 00 
square degrees in its 4 survey regions (e.g. IFu et al1l2008l l. 

In the flat-sky approximation, E and B modes are de- 
fined as 



(21) 



and inversely related as 



that an imperfect survey geometry causes mode-mixing or 



71 (n) ± n2(n) = j ^ [Et' ± ^Bt'] e^(^-^'^^\ (22) 

where the vector n is a fiat-space two-dimensional vector 
that approximates the three-dimensional vector n in Eq. fi]) 
around some reference point; e.g., if the coordinate origin in 
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flat space is taken as the north pole, the position vector of an 
arbitrary point in the vicinity of the north pole is specified 
by two-dimensional vector n = ^(cos sin ip) = {6x,0y)- 
The vector k is the corresponding wavenumber in the flat- 
space coordinate, and (/Sk is defined as k = fc(cos ipk, sin Lpi^). 

Like in the full-sky formalism, we again define pseudo 
E and B modes in a flat-sky limit as 



dn„ [W^(n)(7i(n)±i72(n))]e 



(23) 



where W represents an arbitrary mask field. We take into 
account the mask effect in two steps: a square boundary 
covering observed survey regions and the mask inside the 
square. Even for a complex survey geometry, such a square 
shape geometry can be obtained by using the zero-padding 
method such that the square-shape region encloses the whole 
region of data. 

The pseudo £/_B-mode coeflicient S^^' ± i_B^^* for a 
square field of a side length L relate to the underlying true 
coefficients E^""^ ± iS^^* as (|Bunnll2002l ) 



TTifiat I ■ oflat 



(S^f* ± iB^?')5'k_k'e*^'^''' , (24) 



where Sk is the Fourier transform of the survey window 
function defined as 



^k = 



sin(fc^L/2) sm{kyL/2) 



kx /2 



ky/2 



(25) 



and kx(ky) axe a::(y)-components of k. The window function 
Sk approaches as k goes to 0. 

The shear power spectra in fiat-sky approximation is 
defined as 

{El'^'-E^'*) = (27r)2 5|,(k + k')Cf^ 



/ T-ifiat T-ifiat'i=\ 

(iJk i'k' ) 

riflat oflat* 



(26) 



Here (5|)(k) is the two-dimensional Dirac delta function. 
From Eq. (|24p the pseudo power spectra for the square area 
relate to the true one as 



{2^Y 

Ck'' = L"'(|Bk"T> 



/[cos^(2v9k'k)C^^ +sin^(2(/?k'k)C;^" 



{2-kY 



5k-k'[sm'(2v9k'k)C|;^ + cos2(2v9k'k)C,^- 



(2^) 



r5k-k'[cos^(2vPk'k) - sin^(2v5k'k)]C|;' 



(27) 



the equation (|27p can be approximated by the following al- 
gebraical relations via the mode-mixing matrix: 



5 ^flat 
kk'^-^fc' ■ 



(29) 



Non-zero components of the mode-mixing matrix are given 
using the power spectrum of the power spectrum of survey 
window 5k (eq-[2H]) a-s 

J. _f SEE, EE _ j.,fSBB,BB ^ -^(cos2) 
^^^kk' - ^'^kk' - '^kk' , 



Ti/fS EE , BB 
-''-'kk' 

j.rSEB,EB 

where 



-^{cos2) 
"-"kk' 



^rSBB,EE 
-''-'kk' 



-(sin2) 
•""kk' ) 



-^{cos2) -^(sin2) 
'-'kk' ~ '-'kk' 1 



(30) 



-^{sin2) 
i^kk' 



' IT 27r 



''«+"''^dfc" 2. 

— 5k_k" cos (2v9k"k), 

fe;, -tt/l 



-5k-k" sin (2<^k"k) 



(31) 

Note M-^^^'^^ = = M^BB^EB ^j^g full-sky approach. 
Since the discrete Fourier decomposition is limited by the 
resolution of survey size L, we approximate the mode-mixing 
matrix as an average of the survey window function multi- 
plied with either of cos^ (2(^kk' ) or sin'^ (2(/3kk' ) over a single 
pixel of area (27r/L)^ around the vector k' in Fourier space. 
The mode-mixing due to the survey window is important 
particularly for large-scale modes {k '--^ 2n/L). 

We further need to include the masking effect inside the 
square boundary. As in the full-sky approach, the Fourier 
coefficients of the shear field with the inside mask K{n) is 
related to those without the mask as 



T-ifiat 1 ■ r^flat 



«^*±iB«f') ±2lVk'k, 



(32) 



where the convolution kernel ±2Wk'k is given in terms of the 
Fourier transform of the mask function, A'k, as 



2W^k 



(33) 



In the equation (|32} we use the discrete summation, rather 
than the infinite-range Fourier transform. The summation 
runs over the Fourier-space grids where the resolution of 
each grid is specified by the square size L as {2-k/L)'^. The 
number of grids, N^^^^, needs to be specified by an observer 
to achieve the desired coverage of multipole range. 

The pseudo spectra C^^* for the masked field inside the 
square boundary is 



flat 
k' I 



(34) 



where <^k'k = V'k' — </^k and 5k is the power spectrum of S'k 
defined as 



(Sk50 = (27r)"4(k + k')5k. 



(28) 



The pseudo spectrum in the equation (|27[1 depend on 
wavevector k rather than its length. We later introduce the 
azimuthal angle average over ipu, in reconstructing the true 
spectra Ck that depend only on the length of wavevector k. 

Although the relations between the pseudo and under- 
lying spectra involve an infinite-range integral, we find that 



The mode-coupling matrix for the masking effect is given by 

KBB,BB r-2i^ 2/r, \ ior\ 

= L /Ck-k' cos (2<^k'k), (35) 



M, 



= M, 



KEE,BB 

kk' ^ ^"kk 

KEB,EB T-'^^ r 2/r,.- \ .,:„2 



KBB,EE T-'itr ■ 2/r, \ ina\ 

= L /Ck-k' sm (2i^k'k), (36) 

M^w"'"^" = -f'"''^k-k'[cos^(2</3k'k) - sin^(2iy3k'k)] (37) 

where A3k is the power spectrum of the mask field inside the 
square boundary 



(38) 
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Combining the equations (|29p and (|34|) . the pseudo spectra 
C^^* is related to the true underlying spectra as 

= ^^^k" Mk"k'C^?*. (39) 

k" k' 

We compute the binned power spectra: 

feet 

Cr^^PbfcC^*, (40) 



where the summation corresponds to an azimuthal angle 
average of k, and runs over all the Fourier modes satisfying 
the condition |k| £ [ki,, fc^+i]. The operator Pi,k is defined to 
average dimensionless power within the bin of scales 



^bk — — , 

and its reciprocal is 
2tt 



Q 



fc2 ■ 



(41) 



(42) 



The quantity vi, is the number of Fourier modes available 
for the 6-th multipole bin, and is approximately given for a 
mode of A: S> 1/1/ as 



Vb = Y^ 27vkik[^+'^ - k 



(b) N 
low / 



(43) 



Therefore, as in the full-sky approach, the underlying 
power spectra can be estimated by deconvolving the pseudo 
power spectra in the flat-sky approximation: 



feefc' 



C^' ~ (M-^),,, Y Pb'kiCi""' - (NOmc), 



(44) 



where the mode mixing matrix is given as 

Mbb' =^M^b"Mg„b'. (45) 

b" 

with 

keb k'eb' 
k k' 

feG6 fe'efc' 

Mfft, = ^ftfe ^ Ml^k'Qfc'i,'. (46) 



3.1 Simulating the Lensed Sky 

We construct two types of simulated shear fields: flat-sky 
non-Gaussian fields from ray-tracing simulations to test the 
flat-sky method and full-sky Gaussian fields to test the full- 
sky deconvolution method. 



3.1.1 Full Simulations in Flat-Sky Limit 

We use the shear m aps constru c ted fr om ray-tracing simu- 
lations developed bv lSato et alj |200^ and test the flat-sky 
deconvolution method. The ray-tracing simulations are per- 
formed by placing N-body simulation boxes with two differ- 
ent sizes (240/i~^Mpc and 480/i~^Mpc at a side length) to 
cover a light cone of angular size 5° x 5° from z=0 to 3.5. The 
N-body simulations a re based on th e parallel Tree-Particle 
Mesh code Gadget-2 (|Springelll20"05l ). The ray-tracing sim- 
ulation follows the standard multiple lens plane algo rithm 
(iJain. Seliak fc Whit3l200(]| :[H amana fc Melliei)[200ll ): sep- 
arating the comoving distance between an observer and a 
source by fixed intervals of 120/i~^Mpc, computing the pro- 
jected density contrast at each lens plane along the line-of- 
sight, we obtain the two-dimensional gravitational potential 
field related to the density contrast via Poisson equation. In 
this paper we use the simulations with a single source red- 
shift z — 1. The initial ray directions are set on 2048 x 2048 
Cartesian grids (resolution 300/2048 ~ 0.15 [arcmin]) and 
the ray positions at each lens plane are traced backward to 
the source plane via the lens equation. We obtain the shear 
field from the Jacobian matrix describing deformation of an 
infinitesimal light ray bundle. It should be noted that the 
obtained shear field does not obey the periodic boundary 
condition due to the projection to a single observer point 
through the different simulation boxes. The fundamental 
mode of the shear field is given by Zf — 27r/L — 72 with 
L = (5/180)7r. We use 1000 realizations of the shear maps 
in order to sufficiently reduce statistical scatters in the mea- 
sured power spectra. 

We adopt the concordance ACDM model: the present- 
day density parameters of matter, cosmological constant, 
and baryon are Qn, = 0.238, = 0.762, and Qb = 0.042, re- 
spectively; the Hubble parameter Ho = 73.2 km s~^ Mpc~^; 
the spectral index Us = 0.958; the rms density fluctuations 
in a sphere of 8/i~^Mpc r adius as = 0.76 use d for the power 
spectrum normalization (|Spergel et al.ll2007l ). 



We use the equation (|44p to reconstruct the power spectra 
from the pseudo spectra in the flat-sky approximation. 



3 SIMULATIONS OF SHEAR MAPS 

We use simulated maps of shear flelds to test the pseudo- 
spectrum technique for reconstructing the shear power spec- 
trum. We address how the method works for the case of 
lensing maps, where the lensing power spectrum has distinct 
shapes from the CMB power spectrum: the lensing spectrum 
has greater amplitudes at smaller scales. We also study the 
validity of the fiat-sky approximation for recovering the in- 
put power spectrum down to very small angular scales in 
the presence of the finite-sky survey and masking effects. 



3.1.2 Gaussian Simulations m Full-Sky Limit 

To test the full-sky formalism (Sec. l2.l| ). we co nstruct Gaus- 
sian shear fields using the HEALPix software (|G6rski et al.l 
l2005h . We assume the survey geometry to be a rectangu- 
lar shape covering 2000 square degree (see the top panel of 
Fig-ttl- The resolution parameter of HEALPix pixelization 
Nside is set to be 1024 corresponding to ~ 4 arcmin for the 
pixel scale and the maximum multipole /max ~ 2000. 

The spectrum of the shear flelds is of _E-mode only as- 
suming the concordance cosmological model same as the 
N-body simulations (Sec. |3TTTT|) . In the weak lensing limit, 
the -B-mode shear power spectrum Cf^ is equivalent to the 
power spectrum of lensing convergence field, the weighted 
mass density field projected between source galaxies and an 
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observer alon e the line-of-s ight. Assuming the Limber's ap- 
proximation ( Limbeij|l954l ). the i7-mode power spectrum is 
simply expressed as a weighted projection of the the matter 
power spectrum Ps{k): 



dr r-^W^^{r)Ps{k^ --z 



(47) 



where r is the comoving angular diameter distance, and we 
throughout this paper assume a flat universe. Note that the 
comoving distance is given as a function of redshift, i.e. the 
distance-redshift relation r = r[z). The lensing efficiency 
function, Wahif), is given by 



WGL(r) 



dr'G{r')'- 



(48) 



where G{r)dr = p{z)dz is the normalized source distribution 
and 

p{z) oc Sd{z - Zs) (49) 
with Zs ~ 1. 

3.2 Simulating Realistic Masks 

We consider three types of masks typical for ground-based 
imaging data analyses: point sources, saturation spikes, and 
bad pixels. First masked regions have circular shapes to 
model "bright stars" that are randomly distributed in the 
shear map. The radii of the bright star masks are randomly 
chosen ranging from 0.2 to 2 arcminutes. When the radius 
of the circle is greater than 1 arcmin, we add a rectangular- 
shape mask to model "saturation spikes" around a bright 
star. The mask has the size of 0.2r x 5r (r is the radius of 
masked circle) with the same center as the circle. The dif- 
ferent masked regions are allowed to be overlapped. We also 
mask all the pixels in a j/-direction row randomly selected 
with 5% probability. About 25% fraction of the area is to- 
tally masked in each simulation realization, which is typical 
for a g round-based im aging survey such as the Subaru Tele- 
scope (|Hamana et al.l [20031 1. The configuration of masked 
regions is demonstrated in the bottom panel of Fig. [1] 

Since the pixel size of the ray-tracing simulations (0.15 
arcmin) is smaller than the bright star mask, the mask func- 
tion K{n) is set as if = when the center of the pixel is 
masked, otherwise 1. The middle panel of Fig. [T] shows a 
typical configuration of the simulated map with mask, while 
the spatial resolution is degraded in the plot. On the other 
hand, the Gaussian simulations have a large pixel size (4 
arcmin) compared to the bright star mask, the mask func- 
tion K{n) has a fractional value. Besides this difference the 
masking configuration is kept same in all the realizations for 
both the full-sky and flat-sky simulations. 



3.3 Intrinsic Ellipticities 

The intrinsic ellipticities of source galaxy shapes generate 
a white noise contamination to the power spectrum mea- 
surement assuming random orientations between different 
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Figure 1. Top panel: A representative realization of the full-sky 
simulations assuming that the shear field is Gaussian. The simu- 
lated survey geometry is a rectangular shape with 2000 sq. deg. 
area. Middle panel: Similarly, but for the flat-sky approximation, 
which is done by ray-tracing simulations based on N-body sim- 
ulations for the concordance ACDM model. The simulated map 
has a square shape with 25 sq. deg. area (the corresponding area 
is displayed by the square line in the upper panel). Bottom panel: 
A typical configuration of masked regions in a patch of 20 X 20 
arcmin^ (the corresponding area is displayed by the square line 
in the middle panels). The masked regions occupy about 25% of 
the total simulated area. 
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where is the rms eUipticities per component and fig is the 
mean number density of source galaxies. In the following we 
assume (Te = 0.22 and fig = 30 arcmin"^. 

After masking the simulated shear fields, we add a 
Gaussian random noise with the variance of cr^ /(ngJlpix) 
to each pixel, where Qpi^ is the pixel size. For Gaussian 
simulations, we add a effective noise with the variance of 
^^/(/pixngfipix) to pixels inside the survey area, where /pix 
denotes the unmasked fraction roughly equal to 0.75. 



4 RESULTS 

4.1 Power Spectrum Reconstruction 

Fig.[2]shows the results of the shear spectrum reconstruction 
in the flat-sky approximation using 1000 realizations of the 
ray-tracing simulation. The left panel shows the noiseless 
results, while the right panel shows the results with intrinsic 
ellipticity. The shot noise spectra are subtracted from the 
measured spectra using the theoretical prediction (Eq. [50])) 
in this plot. In reality the shot noise spectra can be estimated 
by erasing the coherent cosmic shear signals in each galaxy 
shape, e.g. by remeasuring the shear spectra after randomly 
rotating orientations of each galaxy. 

Middle panels in Fig. [2] plot the i5-mode power spec- 
trum. The solid curve denotes the input power spectrum, 
while the dashed curve is the convolved spectrum, i.e. the 
measured power spectrum from the simulated shear maps 
without correcting for the survey geometry effect. The con- 
volved spectrum systematically underestimates the underly- 
ing amplitudes over all the scales. On angular scales larger 
than a typical scale of masked areas (small I), the amplitude 
offset is roughly determined by the square of the masked area 
fraction i.e. 0.75 x 0.75 ~ 0.56. On the other hand, the offset 
on larger I is given by the area fraction, 0.75. The modes, 
which angular scale is larger than the mask size, is fraction- 
ally affected by the mask and thus the amplitude decreases 
proportional to the masked fraction. Their power is thereby 
proportional to the square of the unmasked fraction. On the 
other hand, when the scale of modes is mush smaller than 
the mask size, their amplitude is in masked region or 1 in 
unmasked region. As a result, their square value also has a 
value of or 1 and thus the power is proportional just to 
the unmasked fraction (not the square). 

The triangle symbols show the mean pseudo power spec- 
tra over 1000 realizations of the shear maps and the circle 
symbols are the deconvolved one. The error-bars represent 
the expected la statistical errors of the band powers for a 
2000 square degree survey. Since the ray-tracing simulation 
has the area of 25 square degree, we estimate the errors by 
multiplying the measured scatters among 1000 realizations 
with a factor ^2000/25. We set the binning width wider 
than the fundamental mode If — 2n/L = 72: the first 10 bins 
have linearly equal spacing as Z^j^^ — 0.6i/ -I- 1.2Z/6, while the 
remaining 25 bins have logarithmically equal spacing up to 
Zmax = 10000. Our reconstruction method successfully re- 
covers the input i5-mode spectrum for both cases without 
intrinsic noise (left panels) and with noise (right panels). 
The top panel explicitly shows the ratio of the reconstructed 
_E-mode power spectra to the input spectrum. The accuracy 



of the reconstruction achieves a sub-percent level over the 
most ranges of multipoles. 

The lower panels show the B-mode power spectrum rel- 
ative to the i5-mode spectrum. The mask and finite-sky 
effect generates B-mode due to the mode mixing. In fact, 
there is a significant leakage of the iJ-mode into the B-mode 
(~10% of E-mode power) as shown in triangle symbols. We 
find that the leakage disappears to be less than one percent 
after deconvolution (see the filled circle symbols). 

Fig. [3] shows the results of the full-sky approach 
(Sec. 12. ip using 1000 Monte Carlo simulations of Gaussian 
shear fields. The error-bars represent la statistical disper- 
sion for a survey with 2000 square degree coverage. The 
binning width is set to be wider than /f — y^T^/fsky ~ 8: the 
first 10 bins have linearly equal spacing as l^^j^ = 2 + 10 x 6 
with, while the remaining 20 bins have logarithmically equal 
spacing up to Zmax = 2000. The power at the lowest bin 
(12 < Z < 22) is affected by the unseen modes larger than 
the surveyed region and thus the reconstructed power at 
the lowest bin varies with the chosen minimum I and the 
binning width. When the minimum I is taken to be slightly 
larger than If, we find that the reconstructed power becomes 
almost equal to the input power. The accuracy is limited 
by the statistical errors due to insufficient sampling of such 
large-angle modes. The power spectrum without correcting 
for a finite-sky effect (dashed line) is smaller than the input 
power spectrum (solid line). The amplitude offset is roughly 
determined by the square of unmasked fraction times the 
survey area fraction, 0.75 x 0.75 x 2000 square degree di- 
vided by 47r steradian, that is about 0.028. 

Summarizing the results in Figs. [2] and [3] the pseudo- 
spectrum method recovers the input power spectrum at a 
sub-percent accuracy over a wide range of multipoles down 
to I ~ 10". As seen from the top panels, the ratio between 
the true and reconstructed spectra is close to unity over a 
range of multipoles we study. The scatters of the ratio from 
unity imply possible biases in the E-mode power spectrum 
reconstruction at each multipole bin, and are smaller than 
the statistical errors of power spectrum measurement ex- 
pected for a survey of 2000 sq. deg. coverage. The residual 
_B-mode power spectrum is suppressed below a percent level. 

4.2 Correlated Spectrum Errors due to the Mode 
Mixing 

Finite-sky effect makes the errors of the shear power spectra 
correlate between different bins. In this subsection we study 
this effect using the simulated maps. 

Let's start our discussion using Gaussian fields. The sta- 
tistical errors of the binned power spectrum is given, in the 
absence of shot noise, as 

^^(Gau==) ^ (51) 

where /sky is the sky fraction of survey area and vi, is the 
number of modes available around the 6-th bin. 

The left panel of Fig. [4] compares the Gaussian expec- 
tation (|51|l with the statistical errors in the i5-mode power 
spectrum of 1000 Gaussian shear fields. When the multipole 
bin width is taken as Al = 20, 2.5 times wider than the 
the fundamental mode If = -^/tt/ /sky — 8, the errors for the 
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Figure 2. Shear power spectra averaged over 1000 realizations of the ray-traeing simulations. Left panels show the noise- free results, 
while right panels show the results with the intrinsie shape noise. Note that the theoretical prediction for the shot noise power spectrum 
(see Eq. 1501 ) is subtracted from the measured power spectra. Middle panel: The solid lines show the input convergence power spectrum. 
The triangle symbols are the pseudo i?-mode power spectrum directly measured from the masked simulations (without correcting for the 
survey geometry effect), which significantly underestimate the input power spectrum denoted by the solid curve. The dashed lines show 
the theoretical pseudo spectrum obtained by convolving the input spectrum with the mask. The circle symbols show the reconstructed 
B-mode power spectra correcting for the survey geometry effect based on the method in Sec. 12.21 As explicitly shown in the top panel, 
the reconstructed E-mode spectrum agrees with the input spectrum at a sub-percent level over a wide range of multipoles down to 
i = 10*. The deviation is much smaller than the expected statistical errors for a survey of 2000 square degree area coverage denoted by 
the error-bars. The errors are estimated by scaling the measured scatters of the original simulations, which has 25 sq. degree area, with a 
factor ^^2000/25. Bottom panel: The B-mode power spectra before and after correcting for the survey geometry effects, compared with 
the corresponding i?-mode spectra. Note that the B-mode spectrum can be negative after multiplying the mode-mixing matrix with the 
measured E/B-mode spectra, and the cross symbols denote such negative i?-mode values. The reconstruction method suppresses the 
residual S-mode below a percent level shown by the circle symbols. 



simulation results are greater than the Gaussian expecta- 
tion. The limit of sky area and mask causes correlated errors 
between neighboring multipole bins and increases the errors 
by 10-20% over a range of multipoles. When the bin width is 
enough wide such as Al = 100 ^ 12.5if , the power spectrum 
errors become similar to the Gaussian expectations, and the 
spectra at different bins become nearly independent. 

The right panel of Fig. |4] compares the error of the re- 
constructed shear power relative to the unmasked conver- 
gence power using the ray-tracing simulations, which con- 
tain the non-Gaussian error due to nonlinear clustering in 



structure formation. At a bin width A? = 200, a factor 3 
wider than the fundamental mode If = 72, the power spec- 
trum errors increase by about 5% on small multipoles, while 
the correlated errors are well suppressed at a wider binning 
Al = 1000 ~ 10/f. Note that the Gaussian errors depend on 
the bin width of multipoles, while the non-Gaussian errors 
do not: a wider bin relatively suppresses the Gaussi an er- 
ror contributiom^ studied in lTakada fc JainI (|2009l ) (also 
see ISato et al.ir2009l ). the non-Gaussian error contribution is 
significant at multipoles I ^ 1000 for the ACDM model. 
Since most of useful cosmological information in the 
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Figure 3. Similar to Fig. [2] but for the Gaussian simulations in full-sky approach using 1000 realizations. Again the method developed 
in Sec. 12.11 reconstruct the B-mode power spectrum at a sub-percent level. The systematic error in the reconstruction is much smaller 
than the expected statistical errors for a survey of 2000 sq. deg. coverage, both with (left-panel) and without (right panel) the intrinsic 
shape noise. 



shear power spectrum resides in the modes around i ~ 1000 
in the presence of the shot noise, iTakada fc Jaiiil (|2009l ) 
showed that even a factor 2 increase in the power spectrum 
errors at the multipoles degrades accuracies of cosmological 
parameters only by 10-20%. Therefore, the correlated errors 
due to the survey geometry effect have an insignificant im- 
pact on parameter estimation. 



5 SUMMARY AND DISCUSSION 

We develop a pseudo-spectrum method for reconstructing 
the cosmic shear power spectra from actual lensing data. 
The observed shear field is limited to be a finite patch of sky 
and furthermore roughly 25% of the survey area is masked 
due to bright stars. We apply for the first time the pseudo- 
spectrum technique developed in CMB studies to the lensing 
field and show that our method successfully recovers the 
shear spectra over a wide range of multipoles from 100 to 
10^ in both full- and fiat-sky approaches. 



We test the fiat-sky method using ray-tracing simula- 
tions assuming a square-shaped survey region. The 25% frac- 
tion of the total area is masked by reafistic configurations 
for a ground-based survey: circular shapes for bright stars; 
rectangular shapes for bright star spikes; zero padding in one 
direction for bad pixels. We show that both the full- and flat- 
sky methods reconstruct the input i5-mode power spectrum 
in sub-percent accuracy and the residuals are much smaller 
than the statistical errors of a power spectrum measurement 
for a survey of 2000 square degree coverage. The residual 
B-mode power spectrum from the E/ B-mode mixing due 
to the imperfect correction of survey geometry is also sup- 
pressed below a percent of the i5-mode power spectrum. Our 
method offers a new means of measuring the cosmic shear 
correlations and separating the E/B modes from an actual 
survey data. 

Although the pseudo-spectrum technique is promis- 
ing, our method still yields sub-percent residual of B- 
mode in the reconstructed power spectra. To further sup- 
press the B-mode spectrum, one has to remove "ambigu- 
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Figure 4. Left panel: Error increment of the reconstructed shear power spectrum due to the effects of the finite sky area and the 
mask. The errors are measured from 1000 realizations of the fuII-sky Gaussian simulations and then compared with the Gaussian error 
expectation (Eg. lSll V The survey geometry effect increases the power spectrum errors, but it becomes insignificant at a wider bin width 
AZ = 100, about factor 10 wider than the fundamental mode of the 2000 sq. deg. survey area, Zf ~ 8. Right panel: Similar to the left 
panel but for the full simulations in a flat approximation. The error of the reconstructed power is compared with that of the unmasked 
convergence power of the ray-tracing simulations including the non-Gaussian errors. Again the mask effect increases the errors, which is 
suppressed at a wider bin width. Note that the non-Gaussian errors do not depend on the bin width. 



ous mo de" that is inevit ably generated in a finite patch 
of sky (|Bunn et al.l 120031 ). The ambiguous mode satisfies 
both the i5-mode (rotation-free) and the B-mode condi- 
tions (divergence- free) and thus contaminates E/B-mode 
spectrum reconstructed using the simple pseudo-spectrum 
method tha t we a dopt. To ehminate such contamination, 
iBunn et all (|2003l ') introduces pure E/B modes that is or- 
thogonal to the ambiguous modes. Pure pseudo Ci estima- 
tor and its optim ization techn i que of sky a podization have 
been developed JSmithl l2006l: fsmith fc Zaldarriagal l2007l : 



iGrain. Tristram fc StompoJ |2009| : iKim fc Naselskvl |2010| ). 
This technique can be straightforwardly applied to the lens- 
ing case. 

In this paper we assume that masked regions are un- 
correlated with the cosmological shear field. In reality the 
regions with large shear is preferentially masked: in a re- 
gion of massive clusters, we cannot obtain a fair sample of 
background galaxies in the central region due to the dense 
concentration of member galaxies, where the shearing ef- 
fect on background galaxies are greater. Masking such a 
crowded region may bias the power spectrum measurement. 
This contaminating effect can be estimated by combining 
ray-tracing simulations with halo catalogs in the underlying 
N-body simulations. This is our future project, and will be 
presented elsewhere. 
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